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Recovery of interior eigenvalues from reduced near 

field data 

Evgeny Lakshtanov* Boris Vainberg^ 


Abstract 

We consider inverse obstacle and transmission scattering problems where the 
source of the incident waves is located on a smooth closed surface that is a boundary 
of a domain located outside of the obstacle/inhomogeneity of the media. The domain 
can be arbitrarily small but fixed.The scattered waves are measured on the same 
surface. An effective procedure is suggested for recovery of interior eigenvalues by 
these data. 
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1 Introduction 

Consider bounded, connected domain O C M'^, d = 2,3, with twice differentiable boundary 
do. The scattering problem is stated as either an obstacle scattering problem 

-Au^^ - = 0, xE A = > 0, . . 

mSc + = 0, xEdO, 

or a transmission scattering problem 

— — Xn{x)u^^ = A[1 — x E A = > 0, (2) 

where n{x) = 1, x E n{x) >0, x E O, n is smooth in O up to the boundary dO, 

and has interior limiting values at the boundary that are not equal to one, i.e., n{x) — 1 has 
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a constant sign when x G dO. Here is an incident wave that satisfies the Helmholtz 
eqnation in or Here 5 is a bounded set where the sources are distributed. We 

will specify S later. Solution satisfies the radiation condition: 


Akr 


T 2 



e 


X I I 

—, r = |x| —)• cx). 
r 


( 3 ) 


The Dirichlet boundary condition in ([T]) can be replaced by the Neumann boundary 
condition. The transmission problem concerns scattering by inhomogeneous media. The 
media may include an obstacle Oi inside of the region of inhomogeneity. Then equation 
(1^ holds in with the corresponding boundary condition on dOi, see e.g. [1]. We 

will assume that Oi = 0, but all the results below can be automatically extended to the 
case C>i 7 ^ 0. 

Interior eigenvalues are the Dirichlet/Neumann negative Laplacian eigenvalues in O 
for the first problem and interior transmission eigenvalues (ITEs) for the second problem. 
The latter are defined as the values of A for which the following system of equations in O 


— Au — Am = 0, 


—Av — \n{x)v = 0, xeO, u,v e 


u — V = 0, X e do, 

f^-S = 0, 


( 4 ) 

( 5 ) 


has a non-trivial solution. Under our assumptions, the set of ITEs is discrete [28] . 

There is a wide literature on the inverse scattering where the scattering signal is 
registered in a specific direction. A classical example of this situation is the backscattering 
problem. The uniqueness of solution of the inverse problem using the backscattering 
data can be found in 0.0 ,[ID],[23], and the recovery of singularities was studied in 
[21] , |3l] , [22] , [26] . In all the papers above, it was assumed that the echo data is available 
for incident waves coming from all the directions. 

There are important applications where an observer has access to the obstacle only 
from one side. Also the incident waves can be often emitted only from a bounded region, 
and not from infinitely remote points as in the classical backscattering problem. Recently, 
such a potential scattering problem (i.e., problem (]2|) with a potential that is smooth in 
has been studied by Rakesh and Uhlmann [25] in a non-stationary setting. They 
assumed that the incident waves are emitted from points x varying in some sphere. The 
authors show uniqueness for spherically symmetrical potentials. Christiansen [1] stud¬ 
ied the backscattering obstacle problem when the scattering amplitude was known in a 
neighbourhood of the reflected ray. The author obtained a uniqueness result in the class 
of strictly convex obstacles. 

In this article we consider the scattering problems ([T|) and ([2|) when the incident waves 
are emitted from a smooth surface S that is a boundary of a bounded domain B located 
outside of O, i.e., R O = 0. We assume that the receivers are also distributed over the 
same surface S, i.e., the following data is available: 


M*”" emitted from S] 


( 6 ) 
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The key element of our approach is Lemma 12.11 which states that restrictions of waves 
emitted from S form a dense set in L 2 {dO) due to the fact that S is closed. 

We will show that data ([6]) for values of A in some interval X G R'^ allows one to 
determine effectively the interior eigenvalues of the scatterer that belong to the interval 
X. As a byproduct, this leads to the uniqueness of the solution of the inverse problem when 
the data is given for all A > 0 and the spectrum dehnes the scatterer uniquely. The latter 
assumption is true for many classes of obstacles (see |32]) and spherically symmetrical 
potentials (see m. 0.0)- We believe that the arguments below allow for the reduction 
of the amount of scattering data in many methods used in inverse problems to identify 
properties of the scatterer, but here our aim is restricted to the recovery of the interior 
eigenvalues. 

Acknowledgments. The authors are thankful to Rakesh, D.Colton and A.Lechleiter 
for useful discussions. 


2 The main result 

From now on, for the sake of simplicity of notations, we assume that d = 3. The results 
and proofs in other dimensions d > 2 are similar. Dehne operators 

C : L2{S)^H\dO), C* : L2{dO) ^ H\S), 

{Cip){x)= - - -ip{y)dSy, {C*y){x)= y - -fi(y)dSy, k = \/A > 0, (7) 

Js p ~ yi Jao p ~ y\ 

where is the Sobolev space of functions on dO or S. We chose s = 1, but could choose 
any s <2. The upper bound is due to the restriction on the smoothness of dO. 

Lemma 2.1. Let X > 0 be not an eigenvalue of the negative Dirichlet Laplacian in either 
of the domains O or B (with the boundary S). Then operators X, C* have dense ranges. If 
do or S is infinitely smooth, then the range of C (C* respectively) is dense in the Sobolev 
space H^{dO) (H^{S), respectively) with arbitrary s G M. 

Proof. Let us prove that the range of C is dense. Obviously, it is enough to show 
that the kernel of the operator C* is trivial. Assume that the opposite is true. Then there 
exists y = y{x), x G dO, such that y ^ 0 and the function 

r ^ik\x-y\ 

u -.= - - -y{y)dSy, x G k = vA > 0, 

Jdo F - y\ 

which is dehned on and coincides with C*y on S, vanishes on S. Since 

(—A — A)u = 0, X ^ do, 

and A is not an eigenvalue of the Dirichlet problem in R, u = 0 on B. Then from the 
equation above it follows that u = 0 on R^ \ (9. Thus u satishes the Helmholtz equation 
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and the homogeneous Dirichlet boundary condition in O. Since A is not an eigenvalue, it 
follows that M = 0 in O, i.e., u = 0 in The latter contradicts the fact that the jump 
of the normal derivative of u on dO is equal to —dvr/i ^ 0. Thus the density of the range 
of the operator C is proved. Similar arguments are valid for C*. 


□ 


Consider the near held operator 

Fs = Fs{X)-. Fs^ = u^^\s, V e L^iS), X ^ D, 

where is the solution of ([I]) or ([2]) with = Cep on dO (with = Cep on O in the 
case of transmission problem (|2])). 

The main result below iTlieorem 12.31) uses the non-physical operator Fs{X) outside of 
a discrete set of values of A. Discrete sets of A are non-essential since the inside/outside 
duality method determines if a point Aq is an ITE based on the behavior of Fs{X) in a 
neighborhood of Aq, not at Aq. The operator Fs{X) is non-physical since Cip represents 
waves coming to S from inhnity, while waves emitted from S have the form C<p. However, 
one can use operator ^^(A) since each function u^^\s = Fs(p can be obtained (and mea¬ 
sured) using approximations by scattered helds on S produced by some waves emitted 
from S. Namely, the following lemma holds. 

Lemma 2.2. Let X > 0 be not an eigenvalue of the negative Dirichlet Laplacian in 
either of the domains O or B. Then for each (p G L 2 {S), one can construct a seguence 
'f’n £ L 2 {S) of the source densities such that Fgip = lim„^o where the limit is 
taken in the space L 2 {S) and is the solution of (If) or ^ with = £'0n- One can 
determine the source densities without a priori knowledge of O if the radius of a ball 
I a; I < 1/e containing O is known. 

Proof. Consider a bounded domain O that contains O and such that dist(i?, O) > 0. 
For example, one can take O = where B^ is the e-extension of B 

and e > 0 is small enough. Without loss of generality, we can assume that the boundary 
of O is infinitely smooth and A is not an eigenvalue of the negative Dirichlet Laplacian in 

C>. 

From Lemma 12.11 it follows that the range of the operator 





is dense in Then the same is true for C. Hence for every (p G L 2 (iS), there 

exists a sequence x/n £ ^2(5) such that Cx/n —t Fp in Fl^O(^QQf Below we consider 
functions Cxp^, Cip, Cxpn-, Cip defined by the corresponding integrals for all x G The 


standard a priory estimate (e.g., [2l]) for the solution u = Cxp^ — Cip of the Helmholtz 
equation in O implies that 


^ C(A)||D'0^ Cip\\jj 3 / 2 i^QQ'^ 10 as 77. y oo. 
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Since O C O, we have that 

II £93 -11^3/2(90)0 as n-)■ cx) 

and 

||£(/9 — £' 0 „||// 2 (c)) —)■ 0 as n ^ 00 . 

The statement of the lemma is an immediate consequence of the last two relations and 
a priory estimates (e.g., [2T]) for the solutions of problems ([I]) and (|2]) (with radiation 
condition at inhnity). 

□ 

Let us dehne the constant a as follows: a = 1 in the case of the Dirichlet boundary 
condition in the obstacle scattering problem and when n{x) < l,a; G dO, in the trans¬ 
mission scattering problem. Let a = — 1 in the case of the Neumann boundary condition 
and when n{x) > l,x G dO. 

For each function (p G L2{S), consider the complex number {aFs(p, (p) and its argument 
in the interval [0,27r). Let us introduce the function 

$(A) = inf £iTg{aFs{X)ip,ip) e [0,271). (8) 

¥56^2(5), 

The following theorem holds: 

Theorem 2.3. Operator Fs{X) given for all X & X C. M"*", except possibly a discrete set 
of points, determines all the interior eigenvalues in the interval X. Namely, Xq &X is an 
interior eigenvalue of an obstacle scattering problem if and only if 

lim inf <h(A) = 0, cr = 1, liminf [27r — <F(A)] = 0, a =—1, (9) 

A^Ao A^Ao 

where A 7 ^ Aq- The same criterion is valid for the interior transmission eigenvalues of 
odd multiplicity (in particular, for simple ones). 

Remarks. 1) Our proof is based on the inside-outside duality principle [5], [27], 
0, 0, nn, m, [T6] and ideas developed in the study of denseness of the far-held 
patterns m- These principles were used recently na.iig.iisi for numeric evaluation of 
the interior eigenvalues for obstacle scattering, inhomogeneous medium scattering and 
electromagnetic scattering. 

2) In the case of the obstacle scattering, one can easily rehne the statement of the 
theorem and hnd the multiplicity of the eigenvalues. This can not be done for the trans¬ 
mission problem where the inside-outside duality principle allows one to observe not all 
the multiple eigenvalues (with a hxed A), but only some of them, and an even number of 
them can be missed (see more details in ng). One also can consider one sided limits in 
([9]) in the case of the obstacle scattering: A —)■ Aq can be replaced by A t Aq if cx = 1, and 
by A 4 , Ao if O' = —1. 

3) While the proof of the theorem uses the lemma below that contains an assumption 
on the location of the spectrum of the negative Dirichlet Laplacian, the theorem is proved 
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without a priori assumptions on this spectrum since criterium ([9]) does not require the 
knowledge of the operator Fs{X) at all the points A. It is assumed there that A G X with 
an arbitrary discrete set being omitted. 

4) We believe that Theorem 12.31 holds when the boundary dO is Lipshitz. We require 
boundary in order to avoid discussions on existence and properties of Dirichlet-to- 
Neumann maps. 

The following lemmas will be needed to prove the theorem above. 

Denote by 

Fi(A),X°“'(A) ; H\dO) L2{dO) 

the Dirichlet-to-Neumann maps for the Helmholtz equation in the interior and the exterior 
of O, respectively. The solutions are assumed to satisfy the radiation condition when 
is dehned. It is well known that these operators are pseudo differential operators of order 
one (see [29]-[30]) and bounded from to s < 2. The restriction on s is not needed 
if do is inhnitely smooth. When A is positive, operator is analytic in A, and Fi has 
poles at eigenvalues of the negative Laplacian in O. The Dirichlet-to-Neumann map 
for the equation (A -|- Xn)u = 0 in (9 is dehned similarly. The normal vector in all the 
cases is assumed to be directed outside of O. 

Lemma 2.4. Operator Fs{X), A > 0, has the following representation: 

Fs = ^C*{F^-F°^*)C ( 10 ) 

An 

in the case of the Dirichlet obstacle scattering, 

Fs = ^£*Xi((Xi)-i - {F°^^)-^)F^C (11) 

dvr 

in the case of the Neumann obstacle scattering, and 

Fs = fc^F, - F°^^){Fn - F°^y\F, - Fn)C ( 12 ) 

47r 

in the case of the transmission scattering problem. 

Remarks. 1) These formulas are direct analogues of the formulas for the scattering 
amplitude in the problem of scattering of the plane waves (see [6] for the obstacle problem 
and [ISl Th.2.3] for the transmission problem). The only difference is that a plane wave 
is dehned by the direction u of the incident wave, and S is replaced by the unit sphere 
5^ = {(u : |ci;| = 1} in that case. The operators £, C* are also slightly diherent in the case 
of the plane waves. We will denote them by £, C*. Then 

C-.L2{S^)^ L2{dO), £yp(x)= [ (13) 

Js^ 

2) Note that operator (X„ —is analytic in A for A > 0, see [16], Lemma 2.1 item 
2. Thus, the right-hand sides in flTU]) - flT^ are meromorphic in A when A > 0. Operator 
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Fs{X) is analytic there since it is defined via the solution of the exterior problem with an 
analytic Dirichlet boundary condition, and the latter problem does not have eigenvalues 
embedded into continuous spectrum. Hence if relations are established for all 

A > 0 except possibly a discrete set, then they can be extended by continuity toq all 
A > 0. 

Proof. Let us prove (HB. Note that = C(p. We will look for in the form of 
the potential = C*fi with an unknown density /x. Since = —C(p on dO, it follows 
that 

pout^sc ^ _poutp^^ p^^sc ^ _p^p^_ 

Since the jump on dO of the normal derivative of the potential £*/i is equal to —dvr/r, 
from the formulas above it follows that fi = — F°^^)Cip. It remains only to put this 

expression into Fs^p = = £*/i, x E S. 

The case of the Neumann condition is treated similarly using the Neumann-to-Dirichlet 
map instead of the Dirichlet-to-Neumann map. The proof of flT^ is similar to the proof 
of the same formula [HI Th.2.3] in the case of plane incident waves. We will not repeat 
it here since it is somewhat cumbersome. □ 

Let F{k) : ^ 2 ( 6 '^) —)■ ^ 2 ( 6 '^) be the far-field operator, i.e., the integral operator whose 
kernel f{9,oj) is the scattering amplitude u^^{k,9) (see ([3])) defined by the plane incident 
wave . 

Lemma 2.5. Let X > 0 be not an eigenvalue of the negative Dirichlet Laplacian in either 
of the domains O or B. Then the following sets coincide 

{aTg{F{k)p,Lp), p e ^ 2 (^ 2 )} = {arg(F 5 (A)(^,(^), p G ^ 2 ( 5 )}, A = 

where F{k) is the far-field operator. 

Proof. We will give a proof only for the Dirichlet obstacle scattering since the other 
cases are treated absolutely similarly. Due to Lemma 12.41 we have 

(F5(^,(^) = ((F-*-Fi)/i,/i), pi = Cp, 

where C is defined by ([7]) . From the remark following Lemma 12.41 it follows that 

{F{k)p,p) = {{F°^^-F,)ft,fr), pi = Cp, 

where C is defined by (IT^ . It remains to note that operators C and C have dense ranges. 

□ 

Proof of the main theorem. Recall that the scattering matrix S{k) is a unitary 
operator in L 2 {S^) of the form 


S{k)=I + 


ik 

dvr 


F{k). 


(14) 
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Its eigenvalues z = Zj{k) belong to the unit circle. Operator F{k) is compact, and its 
eigenvalues converge to the origin. Thus the eigenvalues Zj{k) of S{k) converge to the 
point z = 1 when k is hxed and j —)■ oo. Functions Zj{k) are analytic in G M"*" everywhere 
except the essential point of the spectrum z = 1. 

We will need the inside-outside duality principle [5], [27], [6], [7], [I2], [IE], which 
concerns the behavior of the set of points {zj{k)} when k changes from zero to oo passing 
through different points k^ for which Aq = fcg interior eigenvalues. The principle states 
(omitting some details that are non-essential for our purpose) that 1 ) one half of the unit 
circle (the upper part if a = 1, and the lower part if a = —1) contains at most a finite 
number of the eigenvalues Zj{k) of S{k) for each fixed value of k, and 2) if m > 0 is the 
multiplicity of an interior eigenvalue A = Ag = fcg, then mi eigenvalues Zj{k) of S{k) from 
the half circle indicated above approach the point z = 1 when k f kn, and m 2 eigenvalues 
Zj{k) from the same half circle approach the point z = 1 when k f fcg, where mi + m 2 < m 
(see [m Lemmas 3.1,3.2] for the duality principle in the case of the transmission problem 
with the assumptions on the refractive index imposed in the present paper). 

These statements can be made more precise in the case of obstacle scattering. Namely, 
mi = m, m 2 = 0 in the case of the Dirichlet boundary condition, and mi = 0, m 2 = m 
in the case of the Neumann boundary condition. In general, these rehned statements are 
not true for the ITEs, see [16], but mi + m 2 = m by modulus 2 in the latter case. 

Since operators S{k) and F{k) have the same eigenfunctions with the eigenvalues 
Zj{k), Xj{k), respectively, and 


Zj{k) = l + ^Xjik), (15) 

the inside-outside duality principle can be reformulated in terms of the far held operator 
F{k) instead of the scattering matrix S{k). Namely, from \zj{k)\ = 1 and ffl^ it follows 
that one can replace condition zj —)■ 1 by either argXj{k) 0 if a = 1 (and Zj{k) belong 
to the upper half circle) or by aig Xj{k) 7 vr if a = — 1 (and Zj{k) belong to the lower half 
circle). 

Next, consider the quadratic form {aFip,ip) and choose aTg{aFip,ip) G [ 0 , 27 r). Then 
the convergence of arg Xj{k) can be replaced by the requirement of the existence of hnite- 
dimensional subspaces on which aTg{aFip,ip) converges to 0 when a = 1 or to 27r when 
cr = — 1. For the sake of simplicity, we assume below that cr = 1. The case of cr = —1 is 
treated similarly. Then the inside-outside duality principle implies that A = Ag = fcg is an 
interior eigenvalue of multiplicity m > 0 if and only if there exist mi (m 2 )-dimensional 
subspaces d*! ($ 2 ) on which arg{F(p,(p) 7 0 when k f ko {k f fcg), respectively. The 
transition from Xj fo the quadratic form is obvious: one can choose to be the span of 
the corresponding eigenfunctions of operator F. The inverse transition is based on simple 
geometrical arguments which follow from ffTTl) and unitarity of S{k). These arguments 
can be found in [ISl Lemmas 3.1,3.2]. The paper [TE] is devoted to the transmission 
problem, but these lemmas have a general character and are applicable to the scattering 
by an obstacle as well. It is worth mentioning that the proof of the inside-outside duality 


principle is often based on establishing the corresponding relation for axg{aFip,ip), and 
then the statement of the principle in terms of eigenvalues of S{k) follows from here 

ng.lEi.iia. 

Due to Lemma [2. 5 [ one can replace F in the quadratic form a.Tg{aFip, ip) by Fs- Thus 
Theorem 12.31 is an immediate consequence of the inside-outside duality principle (stated 
in terms of the quadratic form of F) and Lemma [2.51 

□ 
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